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Abstract 

The PMNS matrix displays an obvious symmetry, but not exact. There are several textures pro- 
posed in literature, which possess various symmetry patterns and seem to originate from different 
physics scenarios at high energy scales. To be consistent with the experimental measurement, all 
of the regularities slightly decline, i.e. the symmetry must be broken. Following the schemes given 
in literature, we modify the matrices (9 in total) to gain the real PMNS matrix by perturbative 
rotations. The transformations may provide hints about the underlying physics at high energies 
and the breaking mechanisms which apply during the evolution to the low energy scale, especially 
the results may be useful for the future model builders. 

PACS: 14.60.Pq Neutrino mass and mixing 



I. INTRODUCTION 



The mixing among fermions is one of the most mysterious aspects in particle physics. 
Unlike the mixing matrix for quarks, the mixing among leptons displays an obvious regularity 
which is manifested in the lepton mixing matrix. It is well known that the mixing among 
fermions originates from the fact that the weak eigenstates of fermions (quarks and leptons) 
are not that of the mass Hamiltonian, and the rotation from the weak basis to the mass 
basis results in the mixing matrix pQ. As observed, the structures of the quark and lepton 
mixing matrices are so different, and it implies the mechanisms which determine their mass 
eigenstates would be different. Lam suggests that a higher horizontal symmetry U(l)x SO (3) 
is broken into the tetrahedral A± and nematic Z2 X Z2 x Z2 sectors which correspond to the 
lepton and quark mixing respectively [2] . Definitely, it is only one of the possible structures 
which were discussed in literature. It is believed that there must be a higher symmetry at 
high energy scales and later it is broken during the evolution from high energy to the weak 
energy scale. It is worth pointing out that the Lam's mechanism which determines an A4 
symmetry for the lepton mixing demands # 13 in the mixing matrix to be zero. And most of 
the proposed symmetries would result in the same zero-^13. However, the recent experiments 
of T2K[3], Double-Chooz [HE], the Daya-Bay [SHE] and RENO[9] collaborations all confirm 
that 813 is not zero, but sizable as near 9°. This implies that even though the lepton mixing 
matrix displays an approximate symmetric form, its original symmetry must be broken. 

The most plausible way to break the symmetry is to perturb the matrix to realize a 
practical lepton mixing matrix which is obtained by fitting the data while the unitarity 
of the matrix must be retained. The form of the perturbation may hint us the breaking 
mechanism which is important for understanding the nature. Moreover, in the process of 
perturbing the matrix and comparing with data, we notice that several textures of the matrix 
are disfavored or marginally favored, even though a perturbative rotation would make them 
to be in marginal agreement with data (see the text). That is the breaking mechanism. A 
careful analysis of the breaking (indeed the perturbation) indicates that one may have an 
opportunity to realize what original symmetric texture(s) is more realistic, so would be able 
to trace back to high energy scale physics where the mixing originates. In particular, such 
a study about the patterns of perturbation may be useful for the future model builders. 

As well known, non-zero neutrino masses; neutrino or lepton mixing and relatively small 
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splitting among neutrino masses are the three conditions leading to the quantum mechanical 
phenomena: observable neutrino oscillations [TUl E] . The mixing matrix in the lepton sector 
LfjU v are named as the Pontecorvo [T2"]-Maki-Nakawaga-Sakata [T3] (PMNS) matrix 



f/pMNS = UjU u . (1) 

which is a 3 x 3 unitary matrix and can be parameterized via mixing angles #12, #23, #13 and 
one CP phase S [TO] 



PMNS 



^ Cl2 c 13 S12C13 S\ 3 e lS ^ 

— S12C23 — Ci2S23Si3e i<5 C12C23 — Si2S23Si3e tS S23 c 13 
y S12S23 — Ci2C23Si3e* <5 — C12S23 — S\2C23Si3e lS C23C13 J 



(2) 



where = cos 9ij, = sin 6^-. If neutrinos are Majorana particles, there could be one 
additional matrix diag(e Ql//2 , e Q2//2 , 1), and since it is not revelent to neutrino oscillations 
at all, we ignore it in this work. The mixing angles and Jarlskog invariant Jcp, which 
determines the magnitude of CP violation in neutrino oscillation [TO1 HH [15] are 

T12 = tan 6»i2 = j^J, (3) 

\U e l\ 

T 23 = tan#23 = P4 (4) 

\ U t3\ 

S 13 = sin6> 13 = \U e3 \, (5) 
J CP = lm(U, 3 U: 3 U e2 U; 2 ). (6) 

The recent data indicate that the angle #13 is sizable: 

• KamLAND Global #13 analysis incorporating CHOOZ, atmospheric, and long- 
baseline accelerator experiments indicates sin 2 6y 3 = 0.009lg 007 ( 1 - e - ^13 = 5.444^2 88i° 
) and non-zero 6*13 at 79% C.L. [T6] . 

• T2K At 90% C.L. and for S CP = 0, 4.99°(5.77°) < 9 13 < 15.97°(17.83°) for normal 
(inverted) hierarchy [3]. 

• MINOS With 5 CP = the best fit result is 2 sin 2 6 23 sin 2 6 13 = 0.041±g;gg(0.079±g;g|i) 
for normal (inverted) hierarchy and 9 13 = is disfavored at 89% C.L. [T7] . 

• Double Chooz The early result from Double Chooz reactor electron antineutrino 
disappearance experiment is 3.7° < 9 r3 < 12° at 90% C.L. [4]. The updated results 



TABLE I: The update global fit results of three neutrino oscillation, where Am 2 is defined as 
m 3 — (mf + m|)/2 and 5m 2 = m 2 , — m\. 
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are sin 2 2# 13 = 0.109 ± 0.030(stat) ± 0.025(syst) (i.e. the central value 6 13 = 9.639°) 
and excluding the no-oscillation hypothesis at 99.8% C.L. [5]. 

• DayaBay The Daya Bay collaboration presents the reactor electron antineutrino dis- 
appearance experiment result sin 2 26'i 3 = 0.092 ± 0.016 (stat) ±0.005 (syst) (i.e. the 
central value #i 3 = 8.8°) and non-zero #i 3 with a significance of 5.2 standard deviations 
[6]. Recent updated result is sin 2 2# 13 = 0.089 ± 0.010 (stat) ±0.005 (syst) (i.e. the 
central value #i 3 = 8.7°) with #i 3 = disfavored at 7.7 a [7] [8]. 

• RENO The result from RENO experiment is sin 2 2#i 3 = 0.113 ± 0.013 (stat) ±0.019 
(syst) (i.e. the central value 9 13 = 9.821°) [9]. 

For convenience of discussion, an updated global analysis on neutrino oscillation data 
[18] is re-presented in Table |TJ and we single out the mixing angles and represent them in 
degrees in Table [TTj 

Analyzing the PMNS matrix, one notices an obvious symmetry, but not exact. If writing 
it in an ideal form which has an exact symmetry, there are various textures which have 
different symmetric patterns. In other words, some phenomenologically assigned forms for 
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TABLE II: The mixing angles from neutrino oscillation fit results in |18j (in degree). 



Parameter 


Best fit 


la range 


2a range 


3(7 range 


6*i2 (NH or IH) 


33.6 


32 6-34 8 


31 6-35 8 


30 6-36 8 


013 (NH) 

_L O \ J 


8.93 


8.45-9.39 


7.99-9.80 


7.47-10.2 


013 (IH) 


8.99 


8.51-9.40 


8.01-9.82 


7.51-10.2 


023 (NH) 


38.4 


37.2-39.8 


36.2-42.0 


35.1-53.0 


023 (IH) 


38.8 


37.5-41.0 


36.5-42.0047.5-53.2 


35.4-54.5 


<5(NH) 


194.4 


138.6-244.8 






<5(IH) 


196.2 


149.4-264.6 







the mixing matrix £7pmns explicitly manifest flavor symmetries while the practical form of the 
matrix implies that the symmetric structures should be spontaneously or explicitly broken. 
Synthesizing the proposals for the symmetric textures existing in literature, there are nine 
in total such ansatzes (1) Tri-Bimaximal Mixing (TBM) [19]; (2) Democratic Mixing (DM) 
[20]; (3) Bimaximal Mixing (BM) [21] ; (4) Golden Ratio Mixing-1 (GRM1) [22] ; (5) Golden 
Ratio Mixing-2 (GRM2) [23]; (6) Hexagonal Mixing (HM) [21]; (7) Tetra-Maximal Mixing 
(TMM) [2S]; (8) Toorop-Feruglio-Hagedorn Mixing-1 (TFH1) [2EH2B]; (9) Toorop-Feruglio- 
Hagedorn Mixing-2 (TFH2) [261428] . We list the explicit forms of these patterns in section 
II. 

Some of the matrix forms listed above require zero-#i3 which is in obvious contradiction to 
the newly measured value. It is shown that all those forms can be modified with perturbative 
rotations into the form of real PMNS matrix which is consistent with data. 

In this work, we explicitly show how a perturbative rotation transforms the the matrix 
into a one with a sizable #13 and practical #12, #23- Our numerical analyses are shown via 
several tables and figures. Then we make some discussions in the last section. 
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II. THE SYMMETRIC TEXTURES OF THE MIXING MATRIX 



Here we list all the nine symmetric textures proposed in literature: 



TBM — 



2 
3 

J_ 

^6 
1 

V6 







1 

^3 
1 1 

1 1 

V3 V2 / 



( 



U~DM — 



\ 



BM 



1 
V2 
1 

^6 
1 

V3 

/ J, 

^2 

_ 1 

2 

1 

2 



1 

1 

^6 
1 

^3 
1 

1 

2 



73 / 
\ 







/ 



GRMl 



1 + 



\ 
( 



U GRM2 



\0 
1 

!(i + V5) 

-\ a/5" 



i 

^2 

71 / 



\ 



^5 



1 

V2 



V5 



\ (5 - V5) 

1+^5 







\ 



HM 



2 



[A 



TMM — 



1_ 

2 a/2 
1 

\ 2\/2 
A) 



4\/2 
1+V5 
4\/2 

\ 

1 

x/2 

72 / 



1 

x/2 

71 / 



l-i(l-V^) 
2~72 
l+i(l-v^) 
2~72 



x/2 



-H 1 + 71 



_j_ 

1+ < 1+ 75) 

2\/2 
2i/2 



( 3+\/3 



TFH1 — 



6 

-3+73 
6 

1 

%/3 



1 

^3 
1 

V3 
1 
V3 



6 

3+v^ 
6 

1 

V3 



6 



u 



TFH2 



/ 3+^ 
6 

1 

Vs 

3-^3 



V 



1 \ 

6 

1 1 

1 3+V3 
V3 6 



(15) 



It is noted that our expressions of the symmetric forms listed above may differ from those 
given in literature by a sign or even a phase factor in a row or column of the matrices, but 
obviously, an additional overall phase e ia does not change the physics of the mixing, and 
moreover, our forms is more convenient to be compared with the conventional expression 
Eq.g adopted by the PDG [TO]. 



III. THE MINIMAL MODIFICATIONS TO THESE PATTERNS 



As well known, the eigenstates of weak interaction are not that of the mass Hamiltonian, 
thus for physical processes one should rotate the weak basis into the mass basis. The unitary 
transformation between the two bases is expressed as a 3x3 matrix: the CKM matrix for 
quarks and PMNS matrix for leptons. The PMNS matrix manifests a not-exact regulation. 
It is supposed that the exact symmetric texture is originating from a symmetry at high 
energy scale and breaking it leads to the practical matrix which keeps an approximate 
symmetric pattern. 

Our goal is to break the symmetric matrix by a perturbation. 

Generally speaking, the perturbation can be realized by transforming the symmetric form 
with two different unitary matrices III ■ Vpmns • Ur- The unitary matrices £/l/r are J us t three- 
dimensional rotations and can be a combination of the following matrices which are rotations 
about three independent axes: 
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(16) 



where 8 X , 5 y , 8 Z are arbitrary phases and s x = sinx, c x = cosx and x, y, z are rotation angles. 
Without losing generality, we only consider the minimal modifications. In this scheme we 
let one of Ul and Ur be a unit matrix, and only the another play the role of perturbation. 

In this work, we only carefully analyze the case for the Tri-bimaximal mixing and an 
explicit illustration on the results is presented by tables and figures, whereas the procedure 
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TABLE III: The results of T 12 , T 23 , S 13 and J C p as perturbing TBM. 



of perturbing the rest eight symmetric textures is similar, so we collect corresponding tables 
and figures in the attached Appendices. 

There are 6 possible ways to perturb the symmetric textures: P x ■ Utbm, Py • ^tbm, 
P z • ^tbm, f^TBM " Px, ^tbm • P y , and C/tbm • Pz- We can obtain the real C/pmns by adjusting 
the parameters in P x , P y , P z . In Table III, we show the trigonometric functions of the mixing 
angles, T 12 , T 23 , S13 and the Jarlskog invariant Jcp. 
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IV. NUMERICAL ANALYSES 



In this section, we analyze the numerical results obtained from the formulation derived 
above. In fact, the procedure for perturbing all these nine symmetric mixing patterns are 
analogous, so we take the Tri-Bimaximal mixing as an example and present the correspond- 
ing results of the rest ones in Appendices |Bj 

Our strategy is following: in the equations presented in last section, we let the left side 
Tij be the experimentally measured values which are based on a global fit of the neutrino 
oscillations and listed in Table [TT], while the right side is the formulas we derived by perturbing 
the symmetric forms. Equating the two sides, we obtain several relations between the model 
parameters, meanwhile we take into account the experimental errors. Plotting them in a 
figure (Fig. [Tj for example), we have three curves which respectively satisfy the relations 
for Ti2 P , T^ p , T^ p . With the experimental errors, the three curves expand into three 
contour bands whose boundaries correspond to the error tolerance. We will observe the 
diagrams and see if they have overlapping regions. If there exists a common region(s) for 
the model parameters where all the three equations are satisfied simultaneously, we would 
say, this scheme is plausible, instead, if there is no such a common region, the scheme is not 
successful and must be abandoned. For instance, in the case of P x ■ Utbm, we have 



T ex P = 2 / 1 + sin(2x) cos 5 X 

12 y 5 + 3 cos(2x) - 4 sin(2x) cos 8 X ' 1 ' 

T 23 P = cosx i ( 18 ) 

s lt P = "7| sina; ' ( 19 ) 

where the superscript "exp" refers to the experimental data. Solving these equations, we 
obtain three curves which correspond to relations between the model parameters x and S x 
as shown in Fig. [1} Due to the experimental errors, the curves expand into bands. The rest 
schemes are similar and we will not respectively discuss the results with different perturba- 
tion ansatzes in every detail, but show them in the following sections and appendices. 

For more explicitly demonstrating the fitting effects, we provide the scatter plots. In 
the plots we set #12, #23, #13 and Jcp as horizontal and vertical axes alternatively, then 
mark the experimental data of the corresponding quantities, and each of them spreads into 
a band whose width is 3 standard deviations (3-<r). There is an overlapping region where 
both experimental data are satisfied within 3cxs. Then we plot our theoretical predictions by 



letting the model (perturbation) parameters scan their whole allowed ranges (for example, 
for P x ■ U TBM , 0° < x < 180° and 0° < 5 X < 180°). If the theoretically predicted values 
which are calculated with a given perturbation ansatz (the red dots) fall into the overlapping 
region, it means that the equation about the model parameters has solutions which coincide 
with the data at least within 3-a tolerance. If there are not red dots in the region, the model 
fails to provide a solution, so that does not work at all. Then even though in all the four 
diagrams solutions for the model parameters seem to exist, we have to investigate if the 
solutions provided by the four scatter plots correspond to the same model parameter region. 
Indeed, the answer resides in the curved band diagrams. Whereas, the scatter plots can 
offer some detailed information about the mixing angles and Jcp which will be measured in 
the future experiments. 

A. P x ■ C/tbm 

We present the curved bands of 812, 823, and #13 in Fig. [T] where the perturbation variables 
x and 8 X serve as the two perpendicular coordinate axes. The bands are obtained by fitting 
the data presented in Table [TTJ The red, light red, and pink regions correspond to 1-a, 2-a 
and 3-a ranges of 6*12, and these three regions are divided by red sold, dashed, and dotted 
lines, respectively. Similarly, the blue, light blue, and nattier blue regions are for 1-a, 2-cr, 
and 3-a ranges of 823 whose boundaries are marked by blue sold, dashed, and dotted lines, 
respectively. The 1-a, 2-a and 3-a ranges of #13 are divided by black solid, dashed and 
dotted lines. In this work, all the curved band diagrams are labelled under this convention. 

The scatter plots among the three mixing angles 812, 823, #13 and Jarlskog invariant Jqp 
are shown in Fig. [2J In the perturbation ansatz P x ■ Utbmj Jcp varies in a range ~ 4 x 10~ 2 . 
Many points fall in the 3-a overlapping region of 812 — 823 and #13 — #12 whereas for # 13 — 823 
the points squeeze on a line which is far away from the central value of 823, as long as we 
require the points not to deviate from the central value of #13 by more than 3<rs. It means 
that simultaneously fitting these two mixing angles is difficult with the P x ■ Utbm ansatz, 
at least not very optimistic. 

By the scatter plots, it is noted that 823 does not exceed 45° within the 3-a ranges of 812 
and 6*13. For the P x ■ £/tbm scheme the conclusion 6*23 < 45° does not change in the whole 
perturbation parameter space of x and S x , i.e., x,5 x G (0°,360°). The constraint #23 < 45° 
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x 

FIG. 1: The solutions corresponding to the measured #12, 623, and #13 for P x ■ JTtbm in the parameter 
space of x — b x . The red, light red, and pink regions are corresponding to the 1-a , 2-a and 3-a 
tolerance levels of #12 (data from Table [TT|) which are divided by red sold, dashed and dotted lines 
respectively. The blue, light blue and nattier blue regions are corresponding to 623 are for 1-a, 2-a 
and 3-cr tolerance levels of 623 (Table |TT|) which are divided by blue sold, dashed and dotted lines 
respectively. The 1-a, 2-a and 3-cr ranges of #13 (without a special color mark in the diagram) are 
divided by black solid, dashed and dotted lines, respectively. Similar diagrams for other schemes 
in this work are labeled with these conventions, (including the color convention.) 

can also be seen from a correlation listed in Table II I II as 

tan 2 9 2 3 = 1 - 2 sin 2 13 (20) 

which indicates that zero 6* 13 results in 623 = 45° or vice versa and non-zero #13 requires 623 
to be less than 45°. 

B. P z ■ C/tbm 

The curved bands of 6*12, 623 and #13 in the whole space of the perturbation variables 
z and S z are presented in Fig. [3j It is noted that the curved band for 9 12 overlaps with 
that of 9 13 within 1-a while its overlap with #23 band is 2-crs from its central value with 
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FIG. 2: The scatter plots for #12, 623, O13 and Jcp with the P x • J7tbm anzatz. The central values 
and 3-0" ranges (by fitting the data shown in Table [TTJ) of the three mixing angles are labeled by 
solid lines and dashed lines, and green for horizontal axis and blue for ordinate one, respectively. 
(Color convention.) 

z G (0°,50°) and S z e (0°,360°). This leads to a conclusion that the perturbation ansatz 
Pz ■ ^tbm is more difficult to accommodate the experimental values of three mixing angles 
simultaneously compared to P x ■ ?7tbm- 

In Fig. [4] we present the scatter plots among the three mixing angles 9±2, 623, #13 and 
Jarlskog invariant Jcp- The perturbation ansatz P z ■ £/tbm provides an upper limit 4 x 10~ 2 
for Jcp. There are many points lie in the 3-cr overlapping region of #13 — #12 while for #12 — #23 
and #i3 — #23, our points fall far away from the central value of 6*23. 
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FIG. 3: The solutions corresponding to the measured 6*12, 023, and #13 for P z -Utbm m the parameter 
space of z — 5 Z . Caption is the same as displayed in Fig. [T] 



From T 23 and Si 3 in Table IIIIl we have a correlation 
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which manifests that a zero-#i3 leads to 823 = 45° or vice versa, while a non-zero #13 deter- 
mines 823 > 45°. 



C J7tbm " -fy 

For ?7tbm ' -Py the curved bands of the three mixing angles are shown in Fig. [5j It is noted 
that 0i2, 023 and #13 share an overlapping region within 1-er. It indicates that the perturbation 
ansatz £/tbm • P y provides a plausible scheme to accommodate all the experimental values 
of three mixing angles. 

The scatter plots among the three mixing angles and Jarlskog invariant in the perturba- 
tion ansatz C/tbm • P y are shown in Fig. |6} The 6>i 3 — J CP presents an upper limit of Jqp 
approximately 4 x 10~ 2 . There are large amounts of points lying in the 3-cr overlapping 
regions of 9 12 — 823 and #13 — #23 while for 813 — 812, points squeeze on a line. Even though 
the line deviates from the crossing point of the central values of #13 and 6*i 2 , this line does 
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FIG. 4: The scatter plots for #12, #23) #13 and Jcp with the P 2 • Utbm ansatz. Caption is the same 
as displayed in Fig. [2] 

pass through the \-a overlapping region of #13 — 6x2. 

Whether #23 > 45° or #23 < 45° cannot be determined in this perturbation ansatz and 
the solution points are observed to be symmetric about the horizontal line #23 = 45°. 

In the scatter plot of #12 — #23, our calculations indicate that as 5 y varies in the range 
(0, vr/2) U (3tt/2, 2vr) 6> 23 > 45° whereas S v E (vr/2, 3vr/2) we note 6> 23 < 45°. This relationship 
can be confirmed by scanning the different parameter ranges of S y presented in Fig. [7j With 
the ansatz £/tbm • Py, U e ^, becomes 
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FIG. 5: The solutions corresponding to the measured 612, 023, and #13 for the C/tbm " Py ansatz in 
the space of parameters y — 5 y . Caption is the same as displayed in Fig. [TJ 

In this case, it is easy to get a conclusion that as y G (0, 7r/2), 8 y G (0,27r), and letting 



8 G (0, 2tt), Eq. (22) and Eq. (23) would demand 8 = 8 y . The equivalence between 8 and 8 y 
implies that the CP phase 5 determines whether 6*23 > 45° or 823 < 45° or vice versa. Table 
[Tf] provides 1-a range of S (the best fit value approximately is 190°) located in (n/2,3ii/2), 
thus 823 should be smaller than 45° with the U^bm • Py ansatz. This implication can also be 
derived from the expressions of T 23 and S13 given in Table III as 



2 1 - sin 2 9i3 + 2 a/2 sin 8 13 ^/l-3 sin 2 8 13 cos S y 

tan 823 = — , ; = = • (24) 

1 — sin 9i3 — 2 v 2 sin 13 y 1 — 3 sin 13 cos 5 y 



The relation indicates that 5 y G (0, 7r/2) U (37r/2,27r) requires 6*23 > 45° and ^23 < 45° as 
5 y G (7r/2, 37r/2). Equality 8 = 5 y means that 6*23 > or < 45° can be determined by the CP 
phase 8 or vice versa. 



D. C/tbm • Pz 

Now let us turn to the ansatz £/tbm • Pz- The curved bands of the three mixing angles in 
the perturbation ansatz ?7tbm • Pz is shown in Fig. [8j It is obvious that the curved band for 
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FIG. 6: The scatter plots for #12, #23 , #13 arid Jqp for ?7tbm • Py Caption is the same as displayed 
in Fig. [2} 

#23 overlaps with that of #13 within 1-a, and 612 band shares overlapping regions with #23 
and 6*i3 within two or three as. 

The scatter plots among these three mixing angles and Jarlskog invariant Jqp are shown 
in Fig. [9] Under the perturbation ansatz ?7tbm ■ Pzi the Jarlskog invariant possesses its 
upper limit approximately 4 x 10 -2 . Plenty of points lie in the 3-cr overlapping region of 
#13 — #23 while for 9\2 — #23 and #13 — #12, the points are far away from the central value of 

#12- 

Similarly to the ansatz C/tbm'-Fj/, &z varies in a range (0, 7r/2)U(37r/2, 2ir), #23 < 45° while 
5 Z G (ir/2, 37r/2), #23 > 45°. The whole range ~ 2tt of 5 Z is scanned and we find that if 
z G (0, 1.3) and 5 Z resides in first and fourth quadrants, #23 < 45°; whereas if S z is in second 



and third quadrants #23 > 45° as shown in Fig. 10 If z > 1.3, there exist mirror-symmetric 
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FIG. 7: The solutions corresponding to the measured #12 and 623 for 5 y G (0, tt/2) U (37r/2, 27r) (left 
panel) and S y G (ir/2,3ir/2) (right panel) with y G (0, tt/2) in the perturbation ansatz £/tbm • -Pz- 
The central values and 3-cr ranges (from the fit in Table [il]) of these two mixing angles are labeled 
by solid lines and dashed lines, green for horizontal axis and blue for ordinate one, respectively. 

diagrams to the corresponding ones. In the perturbation ansatz Utbm ■ Pzi U e 3 is 

(Utbm-P z U = ^=e-"'smz, (25) 



and U e3 in the PMNS matrix is shown in Eq. (23). With z G (0,1.3), 8 y G (0,2vr), and 



5 G (0, 27r), Eq. (25) and Eq. (23) lead to 5 = 5 Z . The equivalence between 5 and 5 Z leads 
to that the CP phase S determines whether 823 > 45° or 823 < 45° or vice versa. With the 
fit value of 5 in Table [TT| we observe that 823 is larger than 45° which is different from the 
conclusion made by the Utbm • Py ansatz. This can also be derived from the expressions of 



T23 and S13 given in Table III 



3 — sin 2 6*i3 — 3\/2 sin 6*13-1/1 — | sin 2 6> 13 cos 8 Z 

tan 2 823 = V , , (26) 

3 — sin 2 6>i3 + 3>/2 sin Am /l — § sin 2 6*13 cos^ 



from which, the ranges of S z G (0, 7r/2) U (37r/2, 27t) and 5 2 G (7r/2, 37r/2) determine 823 < 45° 
or > 45°. 



V. DISCUSSION AND CONCLUSION 

The recent experiments determine a non-zero #13 which is in contrary to the prediction 
made by most of the symmetric textures for the lepton mixing matrix. Even though the 
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FIG. 8: The solutions corresponding to the measured #12, 623, and #13 for U^bm-Pz in the parameter 
space of z — 5 Z . Caption is the same as displayed in Fig. [T] 

real PMNS matrix deviates from the symmetric form, an approximate symmetry is obvi- 
ous. Moreover, it is believed that the symmetric texture is resulted in from the physics at 
higher energy scales and is broken during its evolution to lower energy scales. To inves- 
tigate what physics is at high energy scale, we would study what mechanism breaks the 
symmetry. Following the schemes given in literature, we adopt the perturbation to deform 
the symmetric texture into the real PMNS matrix. Because the original symmetry is ap- 
proximately retained, a perturbation may be a suitable choice. In this work we perturb the 
nine given matrix textures which possess symmetric patterns by various ansatzes. Owing 
to the similarity in all the cases, we take the Tri-bimaximal mixing pattern as an example 



to exhibit how this perturbation method applies. We summarize the results in Table [TV 
Various ansatzes by which the three mixing angles receive corrections are carefully analyzed 
and their effects are marked by a tick y/ or a cross x to note if the ansatz is favored or dis- 
favored. The subscripts S, L, SL of 6*23,5, 6 2 zl and 6*235^ 111 Table IV imply that #23 acquires 
values smaller, larger, and smaller or larger than 45°, respectively. And the mark <g> signifies 
the situation in which the three mixing angles acquire corrections and #13 is non-zero, but 
the ansatz cannot make the theoretical values to be consistent with that data in Table |TT1 



18 






013 



FIG. 9: The scatter plots for correlations among 812, O23, #13 and Jcp for C/tbm ■ P z - Caption is 
the same as displayed in Fig. [2j 



We also mark the best perturbation ansatz with a star symbol * in the table. 

Alternative perturbation schemes have also been proposed. Let us still take 
the Tri-bimaximal mixing as an example to illustrate this new scheme: f^TBM — 
J R 2 3(45 o ) J Ri3(0 o ) J Ri 2 (arctan(l/v / 2)) = J R 2 3(45 ) J Ri 2 (arctan(l/v / 2)), inserting perturbation 
matrices between i? 23 and i?i 2 , namely, i? 2 3 ■ P ■ i?i 2 where P is a suitable perturbation 
3x3 matrix. However, such a perturbation ansatz cannot provide feasible mixing angles to 
be consistent with the experimental data. Therefore such schemes are not phenomenologi- 
cally favorable. 

It is observed that the relationship between 6* 23 and 5i (or 5) (i = x, y, z) could fix 
#23 > 45° or #23 < 45°, thus more precise measurements on #23 constrain the range of the 
CP phase. 
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FIG. 10: The scatter plot for correlation between 612 and 623 for 5 Z £ (0,7r/2) U (3ir/2,2ir) (left 
panel) and 5 Z E (7r/2, 37t/2) (right panel) with y E (0, 1.3) in the perturbation mode UtbwPz- The 
central values and 3-cr ranges (from the fit result in Table [il]) of these two mixing angles are labeled 
by solid lines and dashed lines, green for horizontal axis and blue for ordinate one, respectively. 

The equality 6 — 5i (i — x, y, z) is derived as i is constrained in the first quadrant 

1 G (0, 7r/2), but when it is in the second quadrant i G (tt/2, it), the result would not remain 
the same. From our analytical equations of Jcp, Jcp is proportional to sin(2z) sm<5j (for the 
TMM case, this relationship does not exist) and from standard form of the PMNS matrix 
(J2J) and definition of the Jarlskog invariant (|6]), we determine Jqp ~ sin S. We present a 
relation between Jcp(S) and the perturbation parameters in Table |V} 

Our numerical analysis indicates that the Tri-bimaximal model (TBM) is a more favorable 
texture that may accommodate a sizable #13 after a perturbative correction. With the 
perturbation, 623 and 613 deviate from 45° and 0° as required by the data, and it means that 
the \l — t symmetry [29] originally embedded in the neutrino mass matrix is broken by the 
perturbation. Especially the £/tbm • Py provides the most plausible perturbation ansatz for 
the theoretical mixing angles to be consistent with the experimental values in the l-a level. 
This indicates that the most viable correction to TBM is produced by the rotation in the 

2 — 3 plane, i.e. to break the \i — t symmetry by a perturbation. This provides us a clue for 
the model building in the future. 
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Constant Pattern 




PyU 


Pz-U 


U-P x 






TBM 


V (^235) 


X 


V (^23i) 


X 


V ( ^23 Si)* 


V (^23S_l) 


DM 


® 


X 


(8) 


X 






BM 


V (^235) 


X 


V (^23i) 


X 


® 




GRM1 


v (^235) 


X 


V (^23i) 


X 


V (@23SL,) 


V (®23Sl) 


GRM2 


v (^235) 


X 


V (^23i) 


X 


V (®23Sl) 


V (^23S_l) 


HM 


V(^23s) 


X 


\/(^23l) 


X 


(8> 




TMM 


\/(#23s) 


X 




X 






TFH1 


® 


X 


V(^23l) 


X 




V (023 SL) 


TFH2 


V(®23s) 


X 


<8> 


X 




V(023Sl) 



TABLE IV: The list of the results for perturbing the nine symmetric mixing textures with P x , 
P y , and P z . The mark \J denotes that the corrections from the perturbation scheme make all 
the three mixing angles to be consistent with data, the mark x means that in this ansatz, the 
calculated #13 still remains as exact zero, so such ansatz does not work at all, and <g> for the cases 
that the three corrected mixing angles could not be compatible with the data presented in Table 
|IT| simultaneously, #23 > 45°, 6*23 < 45°, and #23 > or < 45° are marked by subscripts S, L, and 
SL, respectively. The best perturbation ansatz is labeled with *. 
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Appendix A: The analytical results of the rest eight constant mixing patterns 



In this section, we list the results of T12, T23, S13 and Jqp after perturbations to DM 
(Table IVT)), BM (Table |VlT|), GRM1 (Table [VIlTl |ix] |x), GRM2 (Table |Xl] |XlI| |Xffll ), HM 



(Table |XIV| ), TMM (Table |XV[ |XVIl |XVH] |XVIH| ) , TFH1 (Table |XIX] |XX[ ), TFH2 (Table 

xlallxxnl). 
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T23 


S13 
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Px-U 


V 


/ 6+3cos(2:r)+3\/3sin(2:r) cos 8 X 
6+3cos(2x)— 3v / 3sin(2x) cos S x 


1 

V2cosx 


1 O 

w g sin x 


^ sin(2x) sin <5 X 


PyU 


1 
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/ 3+cos(2j/)+2V2sin(22/)cos 8 y 
y 3— cos(2j/) — 2\/2sin(22/) cosiSy 








Pz-U 


\ 


j 5+cos(2z) — 2\/6 sin(2z) cos S z 
1 5+cos(2z)+2\/6sin(22:) cos S z 


\/2secz 


73 sinz 


^sin(2z)sin<5 2 


U-P x 


/l+sin(2a;) cos 5 X 
y 1— sin(2a;) cosia; 


V2 








U-Py 


cosy 


1 /5+3cos(2y)+4sin(2y) cos8 y 

2 y 1— sin(2y) cos S y 


^siny 


Y2 sin(2y) cos<5j, 


U-P z 


secz 


1 /5+3 cos(2z)— 4sin(2z) cos<5z 
2\J l+sin(2z) cos S z 




^ sin(22;) sin<5 2 



TABLE VI: Results of T 12 , T 23 , S 13 and J C p after perturbation to DM. 
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1 








U-Py 


cosy 
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U-P z 


sec z 


/3+cos(2z)-2v / 2sm(22;) cos5 z 


72 sinz 


^sin(22:) sin5 2 
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TABLE VII: Results of T± 2 , T23, S13 and Jcp after perturbation to BM. 



Appendix B: The numerical analyses of the rest eight constant mixing patterns 
1. Democratic Mixing 
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TABLE VIII: Results of T 12 after perturbation to GRM1. 
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1 l+sin(2y) cos 8 y 
\J 1— sin(2t/) cos 5 y 
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TABLE IX: Results of T 23 after perturbation to GRM1. 
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TABLE X: Results of 5 i3 and J C p after perturbation to GRM1. 
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TABLE XI: Results of T 12 after perturbation to GRM2. 
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TABLE XII: Results of T 23 after perturbation to GRM2. 
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TABLE XIII: Results of S 13 and J C p after perturbation to GRM2. 
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TABLE XIV: Results of T i2 , T 23 , S n and J C p after perturbation to HM. 
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TABLE XV: Results of T 12 after perturbation to TMM. 
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Y 2 sin(2j/) [(4+2v^) sin <5 !/ + (v / 2-2) cos 5j,] + (2\/2-l) cos(2j/)+2 v / 2+ll 


u ■ 


Pz 


/- sin(2z)(cos 5 z -4sin<5 z )+2v / 2cos(2z)+5 
Y - sin(2z)(4sin<5 2 +cos5 z )+2v / 2cos(2z)+5 



TABLE XVI: Results of T 23 after perturbation to TMM. 



BM 


Sl3 


Px-U 


\ y / '\/2sm(2x) sin S x + 2 (y/2 - l) sin(2x) cos 4 + § (l - 6^/2) cos(2x) + \ (ll - 2>/2) 


PyU 


1(2-^2) 


Pz-U 


^/-2sin(2z)[2(v / 2-l) cos 5 z - v / 2sin5 z ] + (l-672) cos(2z)-2->/2+ll 

4^ 


U-P x 


1(2-^2) 


U-Py 


^^(2-^2) sin(2y) cos^ + ± (l - 2^2) cos(2y) + | (5 - 2^2) 


U-P z 


^ cos S z sin(2z) - 2 ^2 cos(2z) + 3 



TABLE XVII: Results of S 13 after perturbation to TMM. 
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BM 


Jcp 


Px-U 


X r^sin(2a;) (5\/2smS x + 6cos<L) + 2cos(2x)l 


PyU 


4 cos(2y) 


Pz-U 


k \\$m(2z) (5\/2sin(5 2 - 6cos<5 z ) + 2cos(2z)l 


U-P x 


[(3^/2-2) sin(2x) cos 4 + 4cos(2x)] 


U-Py 


^ K 1 + 7%) sin ( 2 ^) cos *w + 2 cos(22/)] 


U-P z 


^ [2\/2 cos <5 2 sin(2z) + cos(2z)] 



TABLE XVIII: Results of J C p after perturbation to TMM. 



BM 


T12 


?23 


Px-U 


1 2[sin(2:r) cos 5 X +1] 
V -sin(2x)cos(5 a; +^3cos(2x)+2 


^/ \ sin(2x) cos(5 x + jV^cos^x) + 1 


PyU 


y/S- 1 


V 


2(14 


V3) sm(2i/)cos5 !/ +v / 3cos(2?/)+v / 3+4 


- 2 ( 1 + V3) sin(2y) cos 5 a - cos(2j/) + ^+4 


Pz-U 


/ 1— sin(2z) cos<5 z 


^2 (2 + V3) 




Y 2(l+-y/3) sin(2z) cos<S z + v / 3cos(2z)+ v / 3+4 


^/ 2(-v/3-l) cos<5 z sm(2z)+^/3cos(2z)~^/3+4: 


U-P x 


/ 2(l+-v/3) sin(2:r) cos<5 I - v / 3cos(2x)+ v / 3+4 




Y -2(1+^3) sin(2x)cos<5 :E +v / 3cos(2a;)+y3+4 


U-Py 


^/2 ( v^-l) sin(2j/) cos 5 y +a/3 cos(2j/) - ^+4 




/2(l+v / 3) sin(2j/)cos<5 !/ +v / 3cos(2?/)+v / 3+4 


^2(2+^) 


7 1— sin(2j/) cos S y 


U-P z 


/ 2 


\ 


/ — sin(2z) cos<5 z +v / 3~cos(2z)+2 


Y sin(2z)cos5 z + v / 3cos(2z)+2 


' %/2[sin(2z) COS <$* + !] 



TABLE XIX: Results of T 12 and T 23 after perturbation to TFH1. 
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BM 


<Vio 




p x -u 


\J e sm(2x)cos4 2V3 ) + 3 


sin(2x) sin 8 X 
6\/3 


PyU 


6 l V<3 


sin(2y) sin 5,, 
12-^/3 


Pz-U 


^/-2( v / 3- 1 ) cos <5 Z sin(2z) - cos(2z) - 


sin(2z) sin S z 


2^ 


12-v/3 


U-P x 


i(V3-3) 


sin(2x) sin 5^ 
12>/3 


U-Py 


•y/- 2(^-1) sin(22/) cos 5 y - V3 cos(2y) - 


sin(2y) sin 5,, 


2^/3 


12-^/3 


U-P z 


y/— sin(2z) cos 8 Z — \/3 cos(2z)+2 


sin(2z) sin <5 Z 




6^ 



TABLE XX: Results of 5 i3 and J C p after perturbation to TFH1. 



Ti 2 



sin(2:r) cos + 1 



-2(1+^) sin(2:r) cos <5^+ v / 3cos(2x)+V3+4 

>/3- 1 



23 



1 y-2( v / 3-i) sin(2:r) cos fc+v / 3cos(2i:)-\/3+4 

^2(2+^3) 
2(1+a/3) sin(2j/) cos <5j / - v / 3cos(2?/)+v / 3+4 
-2(1+a/3) sin(2?/)cos 5 y + v / 3cos(2y)+y3+4 



2[1— sin(2z) cos 5 Z ] 
sin(2z) cos S z +V3 cos(2z)+2 



— sin(2z) cos 5 Z +V3 cos(2z)+2 



2(l+V3) sin(2:r) cos S x - V3 cos(2x)+v / 3+4 
-2(l+V3) sin(2x)cos5 a; +v / 3cos(2a;)+v / 3+4 



V3- 1 



^-2(>/3-l) sin(2j/) cos <5 y +v / 3cos(2j/)-v / 3+4 
^/2(2+V3) 



sin(2?/) cos Sy+1 



- sin(2z) cos S z +V3cos(2z)+2 



-2(l+%/3) sin(2j/) cos<$ B +-\/3cos(2j/)+-\/3+4 

/ 1— sin(2z) cos S z 

V sin(2z)cos5 z +v / 3cos(2z)+2 



TABLE XXI: Results of T 12 and T 23 after perturbation to TFH2. 
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BM 


<J13 




P x -U 


J 2 ( V3-l) sin(2x) cos 5^ - cos(2:r) - \/3+4 


sin(2x) sin S x 




12\/3 


PyU 


i (3 - V3) 


sin(2j/) sin <5 y 


l2\/3 


Pz-U 


•y/ sin(2.z) cos <5 Z — \/3 cos(2z)+2 


sin(2^) sin S z 




6V3 


U-P x 




sin(2x) sin 5 X 


12^ 


U-Py 


J2(V3- 1 ) cos 5 y sin(2y) - cos(2y ) - V3+4 


sin(2j/) sin 8 y 


2y^ 


12^ 


U-P z 


\J sin(2.z) cos 8 Z — \/3 cos(2z)+2 


sin(2z) sin S z 


v/6 


6^3 



TABLE XXII: Results of S13 and Jqp after perturbation to TFH2. 




FIG. 11: The solutions corresponding to the measured 612, #23, and #13 for P x -Udm in the parameter 
space of x — 5 X . Caption is the same as displayed in Fig. [T] 
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270' 



90° 



FIG. 12: The solutions corresponding to the measured #12, &23, and #13 for P z -Udm in the parameter 
space of z — 5 Z . Caption is the same as displayed in Fig. [T] 

b. P z ■ U DM 
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FIG. 13: The solutions corresponding to the measured #12, #23, and #13 for Uj)M'Py in the parameter 
space of y — S y . Caption is the same as displayed in Fig. [T] 

C. Uvm ■ Py 
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FIG. 14: The solutions corresponding to the measured #12, &23, and #13 for Udm-Pz in the parameter 
space of z — 5 Z . Caption is the same as displayed in Fig. [T] 

d. C/dm • Pz 
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FIG. 15: The solutions corresponding to the measured 9u, &23, and #13 for P x -Ubm m the parameter 
space of x — 5 X . Caption is the same as displayed in Fig. [T] 

2. Bimaximal Mixing 

a- Px ■ Urm 
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FIG. 17: The solutions corresponding to the measured #12, O23, and #13 for P z -Ubm in the parameter 
space of z — 5 Z . Caption is the same as displayed in Fig. [T] 

b. P z ■ U B m 
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FIG. 19: The solutions corresponding to the measured 612, 623, and #13 for UsM'Py in the parameter 
space of y — S y . Caption is the same as displayed in Fig. [T] 

c. Ubm ■ Py 
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FIG. 20: The solutions corresponding to the measured #12, O23, and #13 for Ubm'Pz i n the parameter 
space of z — 5 Z . Caption is the same as displayed in Fig. [T] 

d. Ubm ■ Pz 
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FIG. 21: The solutions corresponding to the measured #12, #23, and 6*13 for P x ■ Uqbmi i n the 
parameter space of x — 5 X . Caption is the same as displayed in Fig. [T} 

3. Golden Ratio Mixing 1 

a- Px ■ t^GRMl 
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FIG. 23: The solutions corresponding to the measured 9u, #23, and #13 for P z ■ E/grmi in the 
parameter space of z — 5 Z . Caption is the same as displayed in Fig. [T] 
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FIG. 25: The solutions corresponding to the measured 612, 023, and #13 for J/grmi • Py in the 
parameter space of y — 5 y . Caption is the same as displayed in Fig. [T] 

C- C^GRMl • Py 
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FIG. 27: The solutions corresponding to the measured #12, 623, and #13 for C/grmi ■ Pz m the 
parameter space of z — 5 Z . Caption is the same as displayed in Fig. [T] 

d- t^GRMl • Pz 
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FIG. 29: The solutions corresponding to the measured #12, #23, and 6*13 for P x ■ C/grm2 i* 1 the 
parameter space of x — 5 X . Caption is the same as displayed in Fig. [T} 

4. Golden Ratio Mixing 2 

a- Px ■ ^GRM2 
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FIG. 31: The solutions corresponding to the measured #12, 623, and #13 for P z ■ [/qrm2 in the 
parameter space of z — 5 Z . Caption is the same as displayed in Fig. [T] 

b. P z ■ ^GRM2 
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FIG. 33: The solutions corresponding to the measured 612, 023, and #13 for J7grm2 • Py in the 
parameter space of y — 5 y . Caption is the same as displayed in Fig. [T] 

C- ^GRM2 • Py 



54 




55 




FIG. 35: The solutions corresponding to the measured #12, 623, and #13 for Uqrm2 ■ Pz in the 
parameter space of z — 5 Z . Caption is the same as displayed in Fig. [T] 
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FIG. 37: The solutions corresponding to the measured 812, #23, and 6*13 for P x -Uhm in the parameter 
space of x — 5 X . Caption is the same as displayed in Fig. [T] 

5. Hexagonal Mixing 

a. Px ■ Uiim 
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FIG. 39: The solutions corresponding to the measured #12, #23) and #13 for P z -Uhm m the parameter 
space of z — 5 Z . Caption is the same as displayed in Fig. [T] 



b. P z ■ Ukm 
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FIG. 41: The solutions corresponding to the measured #12, Q23, and #13 for UuM'Py m the parameter 
space of y — 5 y . Caption is the same as displayed in Fig. [T] 

C. UuM ■ Py 
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10° 20° 30° 40° 50° 



z 

FIG. 42: The solutions corresponding to the measured #12, #23) and #13 for Ubm'Pz m the parameter 
space of z — 5 Z . Caption is the same as displayed in Fig. [T] 

d. Ubm ■ Pz 
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10° 20° 30° 40° 50° 

x 

FIG. 43: The solutions corresponding to the measured 9\2, 6*23, and #13 for P x ■ Utmm in the 
parameter space of x — 5 X . Caption is the same as displayed in Fig. [TJ 

6. Tetra-Maximal Mixing 

a. P x ■ [/tmm 
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FIG. 45: The solutions corresponding to the measured 9±2, 623, and #13 for P z ■ C/tmm i n the 
parameter space of z — 5 Z . Caption is the same as displayed in Fig. [T] 

b. P z ■ Utmm 
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10 20 30 40 10 20 30 40 

#13 <?13 



FIG. 46: The scatter plots for 612, 023, #13 and Jcp with the P z -C/tmm ansatz. Caption is the same 
as displayed in Fig. [2j 623 = 45° labeled as horizontal dotted lines in upper-right and lower-right 
panels. 
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FIG. 47: The solutions corresponding to the measured 612, Q23, and #13 for C/tmm • Py in the 
parameter space of y — 5 y . Caption is the same as displayed in Fig. [T] 
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10° 20° 30° 40° 50° 



z 

FIG. 48: The solutions corresponding to the measured 9±2, 623, and #13 for C/tmm " P% i n the 
parameter space of z — 5 Z . Caption is the same as displayed in Fig. [T] 

d- U^mm ■ P z 
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s 

FIG. 49: The solutions corresponding to the measured 612, #23 > and #13 for P x ■ C/tfhi i n the 
parameter space of x — 5 X . Caption is the same as displayed in Fig. [TJ 

7. Toorop-Feruglio-Hagedorn Mixing- 1 

a. P x ■ Utfhi 
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FIG. 50: The solutions corresponding to the measured 612, 023, and #13 for P z ■ C/tfhi in the 
parameter space of z — 5 Z . Caption is the same as displayed in Fig. [T] 

b. P z ■ C/tfhi 
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10° 20° 30° 40° 50° 

FIG. 52: The solutions corresponding to the measured 612, #23, and #13 for J7tfhi ■ P y in the 
parameter space of y — 5 y . Caption is the same as displayed in Fig. [T] 

c. Ur Fm ■ P y 
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FIG. 53: The solutions corresponding to the measured 612, 023, and 6*13 for C/tfhi " Pz m the 
parameter space of z — 5 Z . Caption is the same as displayed in Fig. [T] 

d. C/tfhi • P z 
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FIG. 55: The solutions corresponding to the measured 812, #23, and #13 for P x ■ C/tfh2 i n the 
parameter space of x — 5 X . Caption is the same as displayed in Fig. [TJ 

8. Toorop-Feruglio-Hagedorn Mixing-2 

a. P x ■ £/tfh2 
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FIG. 57: The solutions corresponding to the measured 9\2, #23, and #13 for P z ■ C/tfh2 in the 
parameter space of z — 5 Z . Caption is the same as displayed in Fig. [T] 

b. P z ■ C/tfh2 
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10° 20° 30° 40° 50° 

FIG. 58: The solutions corresponding to the measured 612, #23> and #13 for J7tfh2 ■ P y in the 
parameter space of y — 5 y . Caption is the same as displayed in Fig. [T] 

C U-TFm ■ Py 
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FIG. 59: The solutions corresponding to the measured #12, 623, and #13 for Uqrm2 ■ Pz in the 
parameter space of z — 8 Z . Caption the same as displayed in Fig. [Tj 

d. C/TFH2 • Pz 
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